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Inferring X-ray coronal structures from Zeeman-Doppler 
images 
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ABSTRACT 

We have modelled the X-ray emission from the young rapid rotator AB Doradus (P ro t 
= 0.514 days) using as a basis Zeeman-Doppler maps of the surface magnetic field. 
This allows us to reconcile the apparently conflicting observations of a high X-ray 
emission measure and coronal density with a low rotational modulation in X-rays. 
The technique is to extrapolate the coronal field from the surface maps by assuming 
the field to be potential. We then determine the coronal density for an isothermal 
corona by solving hydrostatic equilibrium along each field line and scaling the plasma 
pressure at the loop footpoints with the magnetic pressure. We set the density to zero 
along those field lines that are open and those where at any point along their length 
the plasma pressure exceeds the magnetic pressure. We then calculate the optically 
thin X-ray emission measure and rotational modulation for models with a range of 
coronal densities. Although the corona can be very extended, much of the emission 
comes from high-latitude regions close to the stellar surface. Since these are always in 
view as the star rotates, there is little rotational modulation. We find that emission 
measures in the observed range of 10 52 8 — 10 53 3 cm -3 can be reproduced with densities 
in the range 10 9 — 10 10 7 cm -3 for coronae at temperatures of 10 6 — 10 7 K. 
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1 INTRODUCTION 

Studies of stellar coronae have been invigorated recently 



et al. 1993; |Schrijver et al. 199% [Brickhouse fc Dupree 1998|) 



by the weahh of results from Chandra and XMM-Nuwlou. 
Some previous obseivaliuiis with RODAT and ASCA had 

suggested that magnetically active stars have coronae whose 



More recent observations of Capella by FUSE, Chan dra and 
XMM-newton have also indicated high densities (Audard 



hotter gas is more extended than the cooler gas (e.g. Giidel 



et al. 1995). Eclipse mapping of binaries also showed stars 
with extended coronal emission, of which the hotter com- 



ponent suffered less rotational modulation (White et al 



(sec 



199C ). This was not, however a universally accepted view 

IcT 



Singh, White fc Drake 1996 



Siarkowski et al. 1996 



ampapa et al. 1996| ) and the physical extent and location of 



the magnetic loops produ cing this emission was very much 
in question (Jeffries 199S). One of the most challenging re- 
sults however came from line ratio studies that suggested 
that coronal densities in the binaries Capella, a Gem and 



44i Boctis are very high, perhaps up to 10 cm (Duprec 
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et al. 2001; Mewe et al. 2001; Young et al. 2001) (see also 



Giidel et al. 2001b| for XMM-Newton results for a range of 
stars). 

These high densities are problematical because, it is be- 
lieved, such plasmas can only be confined in small, com- 
pact, solar-like loops. These loops cannot explain the lack 
of rotational modulation, unless they are located at very 
high latitudes and are always in view as the star rotates 
or as it is eclipsed by a binary companion. Doppler images 
of active s tars do indeed of ten show polar or high-latitude 
spots (see Strassmeier 1996 for a review) , but they also typ- 



ically show spots at all other latitudes. For AB Dor in par- 



ticular, Zeeman-Doppler images are also ava ilable (Donati 
Henry fc Hall 1995| ; |Donati fc Brown 199?| ). These show 



flux at the kilogauss level at all latitudes. Measurements of 
coronal densities for AB Dor are also high, ranging from 
the 3 x 10 10 cm -3 derived from XMM-Newton observations 
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(Giidel et al. 2001b) to the 3 x 10 9 cm 3 in quiescence and 



lz cm J during flaring deriv ed from modelling a flare de- 



10 

cay observed with BeppoSAX (Maggio et al. 2000). This last 
observation is particularly interesting as it follows the flare 
decay over a timescale longer than a rotation cycle. There 
was no rotational eclipse or modulation of the emission, and 
modelling of the flare decay phase suggested that the loop 
structure was small, with a maximum height of only about 
0.3R*. This suggests that the flaring loop or loops must have 
originated at latitudes above 60° (the stellar inclination) in 
order to remain in view. Observations of directed radio emis- 



sion from AB Dor (Lim et al. 1994) also place the emitting 
regions at high latitudes. 

Simply placing the X-ray emitting loops close to the 
pole may not be the whole solution, however. Emission mea- 
sures for these stars are typically very high. Values deter- 
mined for AB Dor at temperatures fro m 10 6 K to above 10 7 K 



( Giidel et al. 2001b) and HST (Vilhuetal. 2001) have upper 



base d on observations with EXOSAT (Collier Cameron et al i 
198d)~lASCA/EUy E flMewe et al. 1996|), XM M-Newton 



and lower limits of 10 



-10° 



when summed over all 



temperature bins. These emission measures require a large 
volume even at the observed high densities. Any model must 
also e xplain any observed ro tational modulation, or the lack 
of it. Kiirster et al. (1997) presented a long-term study of 



the X-ray emission from AB Dor and concluded that a small 
rotational modulation at the level of 5% — 13% was present. 
Coronal models must also be able to explain the presence 
of very massive prominences often observed in active stars. 
In the case of AB Dor, these are typically observed to be 
trapped in the corona at some 3-5 R* from the rotation 



axis ( Collier Cameron fc Robinson 1989a ; Collier Cameron 
Robi rson"l*989b[ ) . In order for these prominences to form 



there must be sufficient material at these heights to sustain 
the radiative instability thought to trigger their formation. 

The aim of this paper is to use the Zeeman-Doppler 
maps of the surface magnetic field of AB Dor to extrapolate 
the coronal field and hence, by assuming hydrostatic equilib- 
rium, to determine the coronal density and X-ray emission. 
By varying the base density and calculating the emission 
measure and the rotational modulation, we can determine 
if it is possible to reconcile the observations of high den- 
sities and high emission measures with the low rotational 
modulation. 



2 FIELD EXTRAPOLATION 

The meth od of extrapolating the coronal field has been de 



scribed in 



Jardine, Collier Cameron fc Donati (2002) 



and 



will not be repeated in de tail here. We use the source sur- 



face method pioneered by Altschuler fc Ncwkirk, Jr. (1969) 



and a code originally d eveloped by van Ballegooijen, Car- 



tledge & Priest (1998). Briefly, we write the magnetic field 
B in terms of a flux function 9 such that B = — VlP and the 
condition that the field is potential (V x B = 0) is satisfied 
automatically. The condition that the field is divergence-free 
then reduces to Laplace's equation V 2 ^ = with solution 
in spherical co-ordinates (r, 9, <f)) 
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Figure 1. A map of the surface radial magnetic field of AB Dor. 
White represents -800G and black represents 800G. Since AB Dor 
is inclined at 60° to the observer, Zeeman-Doppler images pro- 
vide only limited information in the lower hemisphere. In order 
to compensate for this, we have generated this combined surface 
map, with the 1995 map in the upper hemisphere and the 1996 
map in the lower hemisphere. The observable hemisphere of AB 
Dor typically shows one longitude range of predominantly posi- 
tive polarity, while the remaining range of longitudes is of pre- 
dominantly negative polarity. We have chosen the alignment of 
the 1995 and 1996 maps such that they are antisymmetric about 
the equator, with the predominantly positive polarity regions 180 
degrees apart in longitude in the two hemispheres. 



where the associated Legendre functions are denoted by P; m . 
The coefficients ai m and bi m are determined by imposing the 
radial field at the surface from the Zeeman-Doppler maps 
and by assuming that at some height R s above the surface 
the field becomes radial and hence Be(R s ) ~ 0. This sec- 
ond condition models the effect of the plasma pressure in 
the corona pulling open field lines to form a stellar wind. 
Since large slingshot prominences are observed on AB Dor 
mainly around the co-rotation radius which lies at 2.7i?* 
from the rotation axis, we know that a significant fraction 
of the corona is closed out to those heights and so we set 
the value of R s to 3.4R+. 

We find that much of the corona is filled with open field 
lines, and will therefore be dark in X-rays. These open field 
lines originate in two opposite-polarity, mid-latitude regions 
located 180° apart in longitude. The closed field volume be- 
tween these two "coronal holes" forms a torus that extends 
almost over both poles and whose axis connects the two 
"open field" longitudes. 



3 CALCULATING THE DENSITY 

In order to relate our model to the observations, we deter- 
mine the X-ray emission from the closed field regions. As a 
first step, we calculate the pressure structure of the corona 
assuming it to be isothermal and in hydrostatic equilibrium. 
Hence for a stellar rotation rate uj, the pressure at any point 

is p — poe-f 9sds where g s = (g.B)/|B| is the component of 
gravity (allowing for rotation) along the field and 



9(r,0) 



-GMi,/r + uj rsin 9, u) rs\xi9cost 



(2) 



We note that for AB Dor, w = 1.4 X 10 s , a factor of 
50 times greater than the corresponding solar value. At the 
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Figure 2. Emission measure JnldV for an isothermal corona 
at a temperature of 10 D K for solar maximum (solid line) and so- 
lar minimum (dotted line) as a function of the emission-measure 
weighted density. In each case the branching of the curves at high 
densities shows the amount by which the emission is reduced if 
field lines where the plasma pressure exceeds the magnetic pres- 
sure are assumed to have been forced open. 



loop footpoints we scale the plasma pressure po to the mag- 
netic pressure such that po(6,(f>) = RB (6,<fi) where J? is a 
constant. The plasma pressure within any volume element 
of the corona is set to zero if the field line through that vol- 
ume element is open. In order to mimic the effect of a high 
gas pressure forcing closed field lines to open up, we also 
investigate the effect of setting the plasma pressure to zero 
if the plasma pressure is greater than the magnetic pressure 
i.e. where /3 > 1. From the pressure, we calculate the density 
assuming an ideal gas and determine the morphology of the 
optically thin X-ray emission by integrating along lines of 
sight through the corona. 

This technique has been used for some time now to ex- 
trapolate the coronal field of the Sun. Comparisons of the 
result ing morphology o f the corona with LASCO CI im- 
ages (Wang et al. 1997) show good agreement for much of 
the corona, except near the polar hole boundaries where 
the potential field approximation is likely to break d own at 



the int erface between open and closed field regions. Wang 



et al. (1997) found the best agreement with observations 
was obtained with a scaling of po oc -Bo' 9 . The question of 
the optimal scaling is one that is often add ressed in the 
cont ext of the heating of the solar corona (see Aschwanden 



2001 and references therein) and indeed recent work sug- 



gests that scaling laws developed on the basis of detailed 
solar observ ations may be extrapolated to more rapidly ro- 
tating stars (Schrijver & Aschwanden 2002). Given that our 
intent is not to discriminate between heating models, but 
rather to demonstrate that the observed surface field distri- 
butions are consistent with the X-ray observations, we opt 
for the simplest prescription that fits the data. Our choice of 
po cx Bq ensures that regions of high field strength will emit 
strongly. This naturally concentrates much of the emission 
at high latitudes. 

In order to test the model, we extrapolate the field 
from two solar magnetograms taken from the NSO/KP data 
archive, one from solar maximum (1992 Feb 1) and one from 
solar minimum (1996 April 8). With a source surface set to 



1.40 - 




0.4 0.6 
Rotational Phase 



Figure 3. Rotational modulation of the emission measure for 
three different surface maps. The dashed curve is for the 1995 
dataset. The other two curves are for a mixed dataset with the 
1995 map in the upper hemisphere and the 1996 dataset in the 
lower hemisphere. In each case the relative orientation of the two 
maps is different. The dotted line is the case where the regions 
of predominantly the same polarity in each hemisphere are 180° 
apart in longitude; the solid line is the case shown in Fig. |l| where 
they are at the same longitude. 
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Figure_4. Emission measure for the combined surface map shown 
in Fig. hi with a dipole component of different strengths added. 



2.5 R© and a uniform temperature of 10 6 K we determine 
the total emission measure (/ n 2 e dV) for a range of values 
of the scaling constant R. Although the density at the base 
of the corona is the parameter that we can vary with this 
model, this is not necessarily the density that is observed. 
We therefore determine the emission-measure weighted den- 
sity (J n^dV/ J n^dV) and use this as a measure of the coro- 
nal density. Fig. ^ shows the resulting emission measure for 
the two magnetograms. In each case, emission measures were 
calculated both with and without a cutoff for high pres- 
sures. This causes the branching of the curves at high den- 
sities, with the lower branch being the case where a cutoff 
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was imposed. We do not expect this simple one-temperature 
model to reproduce the detailed changes in the magnitude 
and temperature distribution of the Sun's X-ray emission 
through its cycle. Nonetheless, the overall morphology of 
the emission and the observed range of emission measures 
of 10 49 - 10 50 cm 3 are easily repr oduced and corres pond to 



densities of the order of 10 cm (Peres et al. 2000) 



4 MODELLING THE CORONA OF AB DOR 

In the case of the Sun, resolved X-ray observations can be 
used to test field extrapolation models. For other stars how- 
ever we generally have only the emission measure and the 
rotational modulation to give us insight into the field struc- 
ture. In this section we determine the effect of the different 
assumptions of our model on these two quantities. 



4.1 The effect of hidden flux 

Flux can be undetected by the Zeeman-Doppler method if it 
is located in the unobservable hemisphere of the star, or in 
regions of the stellar surface that are dark, such a s the po- 
lar cap. Jardine, Collier Cameron & Donati (2002) explored 
the effect of this hidden flux by adding to the maps of the 
observable hemisphere either a polar field (in the form of a 
global dipole) or the surface map for another year added in 
to the hidden hemisphere (see Fig. |l|). In all these cases the 
structure of the field (and in particular the positions of the 
open field regions) is quite different. These differences are, 
however, unlikely to be apparent in the rotational modula- 
tion and magnitude of the X-ray emission measure which is 
the main observational signature of this field. 

This is demonstrated in Fig. |^ which shows this ro- 
tational modulation for three different surface maps. The 
lowest curve comes from a surface map that has data only 
in the observable hemisphere. There is very little rotational 
modulation, since flux in the lower hemisphere which might 
be eclipsed as the star rotates is missing. The upper two 
curves show the effect of having flux in both hemispheres. 



In one case, the predominantly positive polarity regions in 
the two hemispheres are aligned, and in the other case they 
are 180° apart in longitude (this is the case shown in Fig. 
0). There is little difference in the rotational modulation of 
these two cases. 

Jardine, Collier Cameron & Donati (2002) also exam- 
ined the change in field structure produced by a "pseudo- 
dipole" field hidden in the dark polar cap. This field is not 
a true dipole since the longest field lines (those that emerge 
from the highest latitudes on the star) extend to the source 
surface and so are forced to be open. As the flux in the dipo- 
lar component is increased, it has a progressively greater 
influence on the overall field structure until, once it dom- 
inates completely, the large-scale field is indistinguishable 
from the pseudo-dipole. Beyond this point, increasing the 
polar field strength has little effect on the field structure, 
but it would increase the emission measure. Since for an 
isothermal corona, the density scales as the plasma pressure 
(and hence in our model as the field strength squared) we 
would expect the emission measure to scale as the polar field 
strength to the power four. Fig. ^| shows the variation of the 
emission measure with the strength of a polar field added 
in to the combined surface map shown in Fig. |l|. Adding 
a weak polar field (< 200G) reduces the emission measure 
as it forces some of the high-latitude field that was origi- 
nally closed to become open. Beyond this field strength, the 
addition of the dipolar component increases the emission 
measure. Even with a polar field strength of a few kilogauss, 
however, the emission measure is increased by only a factor 
of 10. Dipolar fields stronger than this would be observable 
in the Zeeman maps. 

Clearly, while the effect of the hidden flux is important 
in influencing the field topology and the location of the open 
field regions from which the wind may escape, the observa- 
tional capabilities we have at present are not sufficient to dis- 
tinguish the nature of this flux. The different cases we have 
examined all show rotational modulations less than around 
20% which is close to the observed range. While adding a 
dipole field could in principal give an observable change to 
the emission measure, the field strengths that can be hid- 
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Figure 6. The effect of the coronal density on the emission measure. Shown are both emission measure images ('vie wed from longitudes 
128° and 38°) and the corresponding rotational modulation. All models have the same magnetic structure as Fig. |s(a)| , a temperature 
of 10 7 K and a cutoff imposed in the emission for (3 > 1. The top panels are for a model with an emission-measure weighted density 
n e = 4 X 10 8 cm — 3 , while the bottom panels have n e = 1.5 X 10 10 cm — 3 . 



den in the dark polar cap are not large enough to make an 
observable difference. 

In the light of this, we keep the same field topology for 
the rest of this paper, using only the combined surface map 
shown in Fig. [[]. 



4.2 The effect of the coronal temperature 

The coronal temperature influences the X-ray emission pri- 
mari ly th rough its effect on the density scale height. As 
=10 6 K the corona is much more con- 



Fig. 5(c) shows, at T= 
fined than at T=10 7 K. Lowering the temperature not only 
lowers the overall emission level (in this case by two orders of 
magnitude) but it increases the rotational modulation. This 
increase (from 5% to 12%) is, however, small because the re- 
gions that emit most strongly are at high enough latitudes 
that they are always in view. In both these cases, the mag- 
netic structure is the same. Hence, for this low-temperature 
model, large loops still exist, but their summit density is 
much lower and hence they are not so bright in X-rays. This 
may of course have implications for the formation of the 
very massive prominences on AB Dor that are observed at 
some 3-5R* from the rotation axis. If the density at this 



height is too low, it may simply not be possible to sustain 
the radiative instability required to form the prominences. 



4.3 The effect of the coronal density 

This is perhaps the most critical parameter to examine as 
it is currently the subject of such debate. We have chosen 
a sample of base densities and calculated the total emission 
measure and the rotational modulation. In addition, we have 
explored the effect of modelling the breaking open of field 
lines when the gas pressure exceeds the magnetic pressure. 
In these cases, we imposed a cutoff in the emission (i.e. the 
density was set to zero) on those field lines where at any 
point along their length the plasma /3, which is the ratio of 
plasma to magnetic pressure, was greater than one. This cut- 
off has little effect in models where the coronal temperature 
is low (T=10 6 K) since in this case the plasma pressure is also 
low. For the higher temperature models, however (T=10 7 K) 
it makes a more significant difference. 

Fig. ^ shows emission measure images of the X-ray 
corona at two different densities. In both cases, there is a 
cutoff in the emission for those fieldlines with f3 > 1. Also 
shown is the rotational modulation of the emission measure. 
These images show clearly that as the coronal density rises, 
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the accompanying rise in plasma pressure forces open more 
field lines and so the volume of the emitting corona shrinks, 
even although the total emission measure is rising. This is 
accompanied by a slight increase in the magnitude of the 
rotational modulation, although it is still well within the ob- 
served range. The effect of imposing a cutoff in the emission 
where the plasma pressure exce eds t he magn etic pressure 
can be seen by comparing Figs. 5(b) and |3(a)|. For both of 



these models, the viewing angle, bas e den sity and coronal 
temperature are the same, but in Fig. |(b)] there is no cutoff 
imposed. As a result, the emitting volume is independent of 
the density and pressure and the emission measure is larger. 
The rotational modulation is small in both cases. 

In order to quantify this behaviour we have run a series 
of models of increasing base density and calculated the ro- 
tational modulation and emission measure. Fig. |^ shows the 
rotational modulation for a family of models at T — 10 7 K 
where a cutoff in the emission is imposed for ft > 1. The 
modulation is remarkably constant for a very wide range 
of coronal densities. It only begins to increase significantly 
at densities in excess of 10 12 cm -3 . At these densities, the 
plasma pressure at the loop footpoints is so large that even 
very low-lying loops are forced open by the plasma pressure. 
As a result, the rotational modulation increases steeply and 
the emission-measure weighted density in fact falls slightly. 
In the case where no cutoff was imposed (or the temperature 
was so low at T = 10 6 K that the cutoff was never needed) 
the rotational modulation is constant at 5% for all models. 

Although the rotational modulation of the emission 
measure varies little with the density of the emitting corona, 
the magnitude of the emission measure increases steeply. 
Fig. ^ shows the results of models at T = 10 7 K (solid) 
and T — 10 6 K (dashed). At each temperature, the curves 
branch into two: the upper branch is the case where there is 
no cutoff, while the lower branch shows how the emission is 
reduced when a cutoff is imposed. In the case where there is 
no cutoff in the emission, the volume of the emitting corona 
is independent of the density and so the emission measure 
essentially scales as n\. This can be seen in the slope of 
the upper branches. As noted in the discussion of the rota- 
tional modulation, an attempt to increase the base density 
much above 10 12 cm -3 actually results in a reduction of the 
emission-measure weighted density, as can be seen in the 
last few data points. 

Clearly, cutting off the emission when the plasma pres- 
sure exceeds the magnetic pressure results in the emission 
measure increasing less steeply with density. As the density 
increases, more and more field lines reach the critical plasma 
pressure and the emitting volume of the corona shrinks. This 
effect can be seen in Fig. ^| where we show the behaviour of 
the density-weighted filling factor / which we define as 



/ 



JnldV 



(3) 



This is the calculated emission measure as a fraction of 
the emission measure of a sphere extending to the source 
surface R s and uniformly filled with plasma at the coronal 
density. For our isothermal model, this filling factor simply 
measures the fraction of the coronal volume that is emit- 
ting. For the case where there is no cutoff in the emission, 
the filling factor is independent of the density and depends 
only on the structure of the magnetic field and particularly 
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Figure 7. Rotational modulation of the emission measure as a 
function of the emission-measure weighted coronal density. 
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Figure 8. Total emission measure as a function of the emission- 
measure weighted coronal density. Models with isothermal coro- 
na.c at 10 6 K (dashed) and 10 7 K (solid) are shown. In each case the 
curves have two branches. The lower branch shows the amount 
by which the emission is reduced if field lines where the plasma 
pressure exceeds the magnetic pressure are assumed to have been 
forced open. This simulates the opening up of field lines to form 
a stellar wind. 



on the location and extent of the open field regions. For 
the models with a cutoff imposed, however, the filling fac- 
tor falls sharply with increasing density as a smaller and 
smaller fraction of the corona contributes to the emission. 
For models calculated at a lower temperature of 10 6 K, the 
pressure scale height is smaller and this filling factor is cor- 
respondingly less ( 2 x 10~ 4 ) and independent of peak den- 
sity. While the behaviour of this filling factor gives a useful 
qualitative understanding of the structure of the corona of 
AB Dor, the numerical values should be treated with cau- 
tion since they refer only to an isothermal plasma. A true 
stellar corona is a multi-temperature plasma, and defining 
a filling factor is a much more complex issue since fine-scale 
s tructure in both temperature and density m ay contribute 
(lAlmlcakv. Brown fe Sweet 19891 ; | Judge 200 ]] ). Nonetheless, 
Klimchuk fc Cargill (2001) show that the commonly-used 
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Figure 9. Density-weighted filling factor as a function of the 
emission-measure weighted coronal density. The corona is as- 
sumed to be isothermal at a temperature of 10 7 K. The lower 
branch shows the reduction in the filling factor if field lines where 
the plasma pressure exceeds the magnetic pressure are assumed 
to have been forced open. 



spectroscopic filling factor (similar to the one defined here) 
is a good measure of the fraction of the volume occupied by 
plasma within a narrow temperature band. 

The observed range of total emissio n measures for AB 



Dor extends from 10 J 



10 b 



(Vilhu ct al. 2001) 



On the basis of the models shown in Fig. H this corresponds 



to densities in the range io a - irr u - 7 cm . Using a model 
with a cutoff imposed for high-pressure fieldlines reduces 
this range to 10 9 ' 4 — 10 10,5 cm -3 . In either case, the im- 
plied densities are high by solar standards, but span the 
value of 10 10 5 cm~ 3 obta ined from XMM-Newton observa- 
tions (Giidel et al. 2001a). In fact a model run for a temper- 



ature of 3 x 10°K appropriate for these observations gives an 
emission measure of 10 53 ' 6 cm -3 for a density of 10 10,5 cm~ 3 . 



5 CONCLUSIONS 

While our calculation of the coronal magnetic field struc- 
ture of AB Dor is based directly on the Zeeman-Doppler 
images, our calculation of the pressure and density struc- 
ture of the corona is more model-dependent. With the very 
simplest of assumptions however (an isothermal plasma in 
hydrostatic equilibrium) we have been able to reproduce the 
observed emission measures, densities and rotational mod- 
ulation in X-rays. The key to our ability to reconcile the 
high emission measures and densities with the low rotational 
modulation lies in the complex nature of the magnetic field. 
The Zeeman-Doppler images show a surface that is densely 
covered in flux. Extrapolating this surface field shows that 
the corona contains loops on all scales and with a range 
of field strengths. In the low-density models, the corona is 
very extended and so shows little rotational modulation. In 
the higher density models, the emitting corona is more com- 
pact, but again shows little rotational modulation since the 
brightest regions are at high latitudes where they are always 
in view as the star rotates. 
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